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Abstract 

We present preliminary results of the analyses of B —* for and B — ► /iET decays (with h = 
71^, K^ 1 ) from a sample of approximately 60 million BB pairs collected by the BABAR detector at 
the PEP-II asymmetric-energy B Factory at SLAC. We find evidence for a signal in B + — ► 7r + 7r , 
and we measure the branching fraction 

B(B + -► vt+tt ) = (4.liJ;J ± 0.8) x 10~ 6 . 

We also measure the following branching ratios and charge asymmetries: B(B + — * K + -k°) = 
(ll.llJ I ± 1.0) x 10™ 6 , B{B + -» tt + K°) = (17.5±i;f ± 1.3) x 10~ 6 , B{B + K+K°) < 1.3 x 10~ 6 
(90% CL), At+ttO = -0.02+°;^ ± 0.10, A K+7r o = 0.00 ± 0.11 ± 0.02, Ar+if° = -0.17 ± 0.10 ± 0.02, 
where the errors are statistical and systematic, respectively. 
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1 Introduction 



The study of B meson decays into charmless hadronic final states plays an important role in the 
understanding of CP violation in the B system. Measurements of the CP-violating asymmetry in 
the 7r + 7r~ decay mode can provide information on the angle a of the Unitarity Triangle. However, 
in contrast to the theoretically clean determination of the angle f3 in B decays to charmonium final 
states |l], § the extraction of a in tt + tt~ decay is complicated by the interference of b — > uW~ 
tree and b — > dg penguin amplitudes. Since these amplitudes have similar magnitude but carry 
different weak phases, additional measurements of the isospin-related decays^, B + — > 7r + 7r° and 
B° — > 7T 7r°, are required to provide a means of measuring a Q. The measurement of the branching 
ratio of the B + — » 7r + 7r° decay is, in fact, a crucial ingredient, since it is a pure tree amplitude 
decay to a very good approximation. Therefore, in this channel direct CP violation, detected as 
a charge asymmetry (.4), is expected to be zero. Moreover, measurements of B — > Kit decays are 
interesting since phenomenological models have been proposed for extracting the weak phase 7 with 
a global fit to the observables ||, || |(|. We also present here an analysis of the B + — ► tt + K° and 
B + — > K + K° decays. The BABAR collaboration has previously published |m] measurements of the 
branching fractions for B mesons decaying into K + ir° and B + — > ir + K°, but no significant signals 
were seen for B + — > 7r + 7r and B + — > K + K° decays. The results reported here are an update of 
these published analyses. 

2 Data Sample 

The data used in these analyses were collected with the BABAR detector at the PEP-II e + e~ storage 
ring during the years 2000 and 2001. The sample corresponds to an integrated luminosity of about 
54 fb -1 accumulated near the T(4S) resonance ( "on-resonance" ) and about 5fb _1 accumulated at a 
center-of-mass (CM) energy about 40MeV below the T(4S') resonance ("off-resonance"), which are 
used for continuum background studies. The on-resonance sample corresponds to (60.2 ±0.7) x 10 6 
BB pairs. The collider is operated with asymmetric beam energies, producing a boost (/?7 = 0.55) 
of the T(4S) along the collision axis. The boost increases the momentum range of two-body B 
decay products from a narrow distribution centered near 2.6 GeV/ c in the CM to a broad distribution 
extending from 1.7 to 4.3GeV/c. 

BABAR is a solenoidal detector optimized for the asymmetric beam configuration at PEP-II and 
is described in detail in Ref. ||. Charged particle (track) momenta are measured in a tracking 
system consisting of a 5-layer, double-sided, silicon vertex tracker and a 40-layer drift chamber 
filled with a gas mixture of helium and isobutane, both operating within a 1.5 T superconducting 
solenoidal magnet. Photon candidates are selected as local maxima of deposited energy in an elec- 
tromagnetic calorimeter (EMC) consisting of 6580 CsI(Tl) crystals arranged in barrel and forward 
endcap subdetectors. In this analysis, tracks are identified as pions or kaons by the Cerenkov angle 
9 C measured by a detector of internally reflected Cherenkov light (DIRC). The DIRC system is 
a unique type of Cherenkov detector that relies on total internal reflection within the radiating 
volumes (quartz bars) to deliver the Cherenkov light outside the tracking and magnetic volumes, 
where the Cherenkov ring is imaged by an array of ~ 11000 photomultiplier tubes. 



Charge conjugate modes are assumed throughout this paper. 
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3 Event Selection, 7r° and K° Reconstruction 

Hadronic events are selected based on track multiplicity and event topology. Backgrounds from 
non-hadronic events are reduced by requiring the ratio of Fox- Wolfram moments, H2/HQ ||, to be 
less than 0.95 and the sphericity [[[(]] of the event to be greater than 0.01. 

Candidate 7r° mesons are reconstructed as pairs of photons with an invariant mass within 3a 
of the nominal 7r° mass [11], where the resolution a is about 8 MeV/c 2 . Photon candidates are 
selected as showers in the EMC that have the expected lateral shape, are not matched to a track, 
and have a minimum energy of 30 MeV. The tt° candidates are then kinematically fitted with their 
mass constrained to the 7r° nominal mass. 

mesons are detected in the mode K° — > K® — » ir + ir and are reconstructed from pairs 
of oppositely charged tracks that form a well-measured vertex and have an invariant mass within 
11.2MeV/c 2 (which corresponds to 3.5a) of the nominal K® mass [p7lf| . The measured proper decay 
time of the K® candidate is required to exceed five times its uncertainty. 



4 B Reconstruction 

B meson candidates are reconstructed by combining a it or a K® candidate with a track h. The 
kinematic constraints provided by the T(4S) initial state and knowledge of the beam energies are 
exploited to efficiently identify B candidates. We define a beam-energy substituted mass toes = 

yl ' E^ — p 2 B , where E^ = (s/2 + • p^)/^, \/s and Ei are the total energies of the e + e~ system 
in the CM and lab frames, respectively, and pi and p# are the momentum vectors in the lab 
frame of the e + e~ system and the B candidate, respectively. An additional kinematic parameter 
AE is defined as the difference between the energy of the B candidate and half the energy of the 
e + e~ system, computed in the CM system. The tti.es resolution is dominated by the beam energy 
spread, while for AE the main contribution comes from the measurement of particle energies in 
the detector. These two variables are therefore substantially uncorrelated. 

However, in the hir° (with h = tt^jK^) final states both the AE and tttes distributions have a 
tail due to imperfect containment of the electromagnetic showers initiated by the 7r°. In this case 
only, in order to reduce this source of correlation and to slightly improve the resolution, we fit the 
B candidate with the energy constrained to the CM beam energy in the two cases of kaon and 
pion mass hypothesis for the track h. For the /wr decay the energy-constrained mass resolution 
is then found to be about 3 MeV/ c 2 from the core Gaussian width of a Crystal Ball[] fit to Monte 
Carlo simulated signal events. For the hK® decay the mES resolution is found to be 2.5 MeV/ c 2 
from a Gaussian fit. For both decay topologies the signal Monte Carlo resolutions are validated 
by comparing data and Monte Carlo resolutions for decays into open charm final states with large 
branching fractions, such as B~ — ► D°p~ , (with p~ — ► 7r~7r° and D° — > K~ir + ) for the hir analysis, 
and B~ -> D°ir~ (D° -> K~tt + ) for the hK° s analysis. 

The AE variable is evaluated assuming the pion mass hypothesis for the track h. Its distribution 
for the signal 7r + 7r° events is described by a Crystal Ball function centered near zero. Since the 
AE distribution has a mean that depends on the track h momentum in the lab frame in the case 
of signal K + tt° events, we also calculate AE with the kaon mass hypothesis (AE(K)) for those 
events. We empirically find that its distribution is described better by a sum of two Gaussians 
with different mean values. For hK® signal events the AE distribution is parametrized as a sum 

2 A core Gaussian with a power law to describe a tail at negative values is called the Crystal Ball function |12| . 
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of two Gaussians centered near zero with the core Gaussian accounting for 95% of the events, 
taking into account the momentum dependance for signal B + — > K + K®. Based on Monte Carlo 
simulated B + — > 7r + 7r° and B + — > TT + Kg events, we estimate the resolution on AE for the core 
Gaussian width to be about 40 MeV and 26 MeV, respectively. Candidates are selected in the 
range 5.2 < th-es < 5.3GeV/c 2 . Different requirements on AE specific to each analysis are then 
applied. 

5 Background Rejection 

The dominant background to these channels is from random combinations of a true ir° (Kg) with 
a track, produced in e + e~ — ► qq continuum events (where q = u, d,s, or c). Another source of 
background originates from B decays into three (or more) light mesons. Detailed Monte Carlo 
simulation, off-resonance, and on-resonance data are used to study backgrounds. For this study we 
select on-resonance data in AE sideband regions defined by the ranges 0.20 < \AE\ < 0.45 GeV for 
hn°, and -0.305 < AE < -0.115 GeV plus 0.075 < AE < 0.265 GeV for hK°. 

In the CM frame the continuum background typically exhibits a two-jet structure, in contrast 
to the isotropic decay of BB pairs produced in T(45) decays. We exploit the topology difference 
between signal and background by making use of two event-shape quantities. 

The first variable is the angle 6 S between the sphericity axes of the B candidate and of the 
remaining tracks and photons in the event. The distribution of | cos 9 S \ in the CM frame is strongly 
peaked near 1 for continuum events and is approximately uniform for BB events. We require 
| cos# s | < 0.8 in the /i7r° analysis, but, given the lower level of background, only | cos# s | < 0.9 in 
the hK® analysis. 

The second quantity is a Fisher discriminant J- |7J constructed from the scalar sum of the CM 
momenta of all tracks and photons (excluding the B candidate decay products) flowing into nine 
concentric cones centered on the thrust axis of the B candidate. Each cone subtends an angle of 
10° and is folded to combine the forward and backward intervals. Monte Carlo samples are used 
to obtain the values of the Fisher coefficients, which are determined by maximizing the statistical 
separation between signal and background events. No requirement is applied on T\ instead the 
distributions for signal and background events are included in a maximum likelihood fit as described 
in the next section. 

On the other hand, B background events tend to peak in mEs> as do signal events, but have 
more negative AE values. They are particularly harmful for the hir° analysis given the poorer AE 
resolution. We use data in the negative AE sideband region to estimate the magnitude of this 
background and Monte Carlo techniques to choose a AE requirement that reduces this background 
to a negligible level. We finally require -0.11 < AE < 0.15 GeV for /i7r° and -0.115 < AE < 
0.075 GeV for hK° s . 

A total of 13661 candidates in the on-resonance data satisfy our hir° selection criteria and 
with the fiK® analysis requirements we select 10668 candidates. These two samples enter into two 
separate maximum likelihood fits. 

The final selection efficiency e is (25.6 ± 1.7)% [(22.5 ± 1.5)%] for B + -> vr+vr [B + -> K + ir°] 
events, while it is (47.5 ± 2.2)% [(47.1 ± 2.2)%] for B+ -» tt+K® [B+ K+K°] events. The 
errors on the efficiencies are statistical and systematic, combined in quadrature. The dominant 
component is due to the imperfect knowledge of 7r° and K® reconstruction efficiencies (5% and 3% 
relative errors, respectively). 
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6 Signal Extraction 



For each topology (hir° and hK®), an unbinned maximum likelihood fit determines the signal and 
background yields ni (i = 1 to M, where M is the total number of signal and background species) 
and charge asymmetries A% = (n^ — nf)/(n^ + nf), where n~ (nf) is the fitted number of i th 
type h~ir° (h + ir°) [h~K® (h + K®)] events. The input variables to the fit are wies, AE, T and the 
Cherenkov angle 9 C of the track from the candidate B decay. The extended likelihood function C 
is defined as 



/ M \ N 

£ = exp I - ^2 TH J II 

V %=\ ) j=i 



M 1 

^ -(1 - qjAi)niVi (xj;di 
,i=l 



(1) 



where qj is the charge of the track h in the j th event. The M probabilities Vi(xj;di) are evaluated 
as the product of probability density functions (PDFs) for each of the independent variables Xj, 
given the set of parameters <3j. Monte Carlo simulation is used to validate the assumption that 
the fit variables are uncorrelated. The exponential factor in the likelihood accounts for Poisson 
fluctuations in the total number of observed events N. 

The parameters for the background ttjes and AE PDFs are determined from events in the 
off-resonance data and in the tties sideband region, respectively. The ttt-es shape is parameterized 
by a threshold function ]l3| /(mEs) oc ttiesVI ~~ x2 ex P[ — £(1 — x<2 )\i where x = mEsA^o and mo 
is the average CM beam energy. The background shape in AE is parameterized as a second-order 
polynomial. The signal distributions have been already described in Sect. ^. 

Events from Monte Carlo simulated signal decays and from on-resonance mEs sideband regions 
are used to parameterize the Fisher discriminant PDF for signal and background events as a 
Gaussian and a sum of two Gaussians, respectively. Alternative parameterizations for obtained 
from off-resonance data (for background) and B~ — > D°ir~ fully reconstructed decays (for signal), 
are used to estimate systematic uncertainties. The 6 C PDFs are derived from kaon and pion tracks 
in the momentum range of interest from a sample of D* + — > D°tt + (D° — > K~7r + ) decays. This 
control sample is used to parameterize the 6 C resolution as a function of track polar angle. 

The results of the fit are summarized in the first column of Table |], where the statistical error 
for each mode corresponds to a 68% confidence level interval and is given by the change in signal 
yield rtj that corresponds to a — 21n£ increase of one unit. We define a signal statistical significance 
as the square root of the change in — 21n£ when the signal yield is fixed to zero. For the 7r + 7r° 
mode, we find a 5.2a statistical significance for the signal. 

In order to increase the relative fraction of signal events of a given type for display purpose 
we choose events passing requirements on likelihood ratios. These likelihood ratios are defined as 
= Yja'Ps/Yli'Pi aim ^-k = Vhj SsPs, where denotes the sum over the probabilities for 
signal hypotheses only, denotes the sum over all the probabilities (signal and background), and 
Vk denotes the probability for signal hypothesis k. These probabilities are constructed from all 
the PDFs except that describing the plotted variable. Figures [j] and || show the distributions in 
mES and AE for events passing all such selection criteria. The likelihood fit projections, scaled 
by the relative efficiencies for the likelihood ratio requirements, are overlaid on each distribution. 
Since the sample projections in toes and AE are obtained with requirements on different likelihood 
ratios, the number of signal events appearing in the two projections are not the same. 
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Figure 1: Distributions of tties and AE for 7r + 7r° events (left) and K + ir° events (right) after 
additional requirements on likelihood ratios, based on all variables except the one being plotted. 
Solid curves represent projections of the complete maximum likelihood fit result; dotted curves 
represent the background contribution. 
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Figure 2: Distributions of tties (left) and AE (right) for ir + Kg events after additional requirements 
on likelihood ratios, based on all variables except the one being plotted. Solid curves represent 
projections of the complete maximum likelihood fit result; dotted curves represent the background 
contribution. 



Table 1: Summary of fitted signal yields, measured branching fraction B and charge asymmetries 
A. The first error is statistical and the second is systematic. For the K + K° mode we quote the 
90% confidence level (CL) upper limits for the signal yield and branching ratio, and give the central 
values in parentheses. 

~Mode Signal Yield g (10"**) A 

vr+vr 11 62+^ ± 11 4.1+^ ± 0.8 -0.02+^ ± 0.10 

K + tt° 149 ±17 ±8 ll.lijjil.O 0.00 ±0.11 ±0.02 

ir+K° 172 ±17 ±9 17.5±i;f db 1.3 -0.17 ± 0.10 ± 0.02 

K + K° < 10 (— 5.6±|;| ± 2.5) < 1.3 (— 0.6±g;? ± 0.3) 



7 Branching Fraction Results 

The branching fractions are defined as 

^) = \ > &> 

£(vr u -» 77) e h7V o ■ N BB - 

1 n hK° 

B{kK ) = B(K° - Kf) ■ B(K° - vr+vr-) e hK% • N bE ' (3) 

where n^o (n hK o) is the signal yield from the fit and e^o (e^o) is the reconstruction efficiency 
for the mode hir° {hK° s ) in the detected tt° (K° s ) decay chain. N BB = (60.2 ± 0.7) x 10 6 is the 
total number of BB pairs in our dataset. B(tt° — » 77), B(K° — > Kg), and B(K® — > ir + 7r~) are 
taken to be equal to 0.98798, 0.5 and 0.6861, respectively pl| . Implicit in the above equations is 
the assumption of equal branching fractions for T(4S) — > B°B° and T(4S) — ► B + B~ . 

Systematic uncertainties on the branching fractions arise primarily from uncertainty on the final 
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selection efficiency and uncertainty on n; L due to imperfect knowledge of the PDF shapes. The latter 
is estimated either by varying the PDF parameters within la of their measured uncertainties or 
by substituting alternative PDFs from independent control samples. In the hit analysis the most 
relevant systematic uncertainties on the signal yields are due to the background ttt-es parametriza- 
tion and Fisher background shape (about 10% each), while for the hK® analysis the AE offset and 
resolution and Fisher signal shape contribute the largest errors (about 4% each). We estimate the 
systematic uncertainty on the signal yields due to the residual presence of B decay backgrounds 
with Monte Carlo techniques and we find that it is negligible compared with the other effects. 

In the case of the ir + ir° final state, we evaluate how the imperfect knowledge of the PDF shapes 
can affect the significance of the signal. We recalculate the square root of the change in — 21n£ 
with n n + w o fixed to zero for the worst case PDF variations and we find a 4.0a statistical significance 
for the signal. 

Systematic uncertainties on the charge asymmetries are evaluated from PDF variations added 
in quadrature with the limit on intrinsic charge bias in the detector (0.01). The small yield of ir + ir° 
channel is the origin of the systematic error on the charge asymmetry (0.10), dominated by the 
PDF variations. 

In conclusion, we find evidence for the decay B + — ► 7r + 7r° and measure a branching fraction of 
B{B+ -► vt+tt ) = (4.r^;£±0.8) x 10~ 6 . We also measure B{B+ -► K + tt°) = (11.1±J;|± 1.0) x 10~ 6 
and B(B + — > ir + K°) = (17.5l^'y ± 1.3) x 10~ 6 , with significant improvements on the errors with 
respect to our previously published results. We do not observe any evidence of direct CP asymmetry 
in these channels, measuring A^+^o = —0.02^026 ± 0.10, Ak+tt° = 0.00 ±0.11 ±0.02, and A^+k = 
— 0.17±0.10±0.02. No evidence of a signal is found for the K + K° final state for which we set a 
90% CL upper limit on the branching ratio of 1.3 x 10~ 6 . 
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